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Abstract 

In this article the error estimation of the moving least squares approximation is provided 
for functions in fractional order Sobolev spaces. The analysis presented in this paper 
extends the previous estimations and explains some unnoticed mathematical details. An 
application to Galerkin method for partial differential equations is also supplied. 
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1. Introduction 


The Moving Least Squares (MLS) approximation was introduced in an early p^er 
by Lancaster and Salkauskas [l| in 1981 with special cases going back to McLain [JQ 
in 1974 and 1976 and to Shepard Q in 1968. For other early studies we can mention the 
work of Earwig [^, 0, 0]- Since, in MLS one writes the value of the unknown function 
in terms of scattered data, it can be used as an approximation to span the trial space 
in meshless (or meshfree) methods. This approximation has found many applications in 
curve fitt ing and numerical solutions of partial differential equations since early nineties 
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The error analysis of MLS approximation was provided by some authors, beginning 
with the work of Earwig 0] which is limited to a univariate case. The connection to 


Backus-Gilbert optimality was studied by Levin 12| in 1998, and later it was used by 




Wendland [la M lla in a more elaborated setting. In Liu et. al. [16| the anal ysis 
is presented for smooth functions in fl Armentano and Duran 171 

proved error estimates in L°° for the function and its first derivatives in one dimensional 
Afterward Armentano 1^ generalized this to multi-dimensional cases but it is 


case. 


still restricted to “convex” domains and Sobolev spaces of order one. One can also find 
an estimation in Han and Meng [i0 | for reproducing kernel particle methods (which is 
related to the MLS approximation) for integer order Sobolev spaces. They assumed a 
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constant bound for the norm of the inverse matrix (matrix A in text) and considered it for 
special cases in one dimension and first order approximations. Note that the role of this 
matrix is very crucial in analysis. The paper of Zuppa 0 is also limited to some specific 
situations. In Wendland 13|, [l^ the analysis presented only for the function in classical 


function spaces. We can also mention the work of Melenk 2l| where the theoretical 
and computational aspects of some meshless approximation methods, including MLS, 
are considered. 

The collocation method based on the MLS approximation is called finite point method. 
An analysis for this method has been presented in 22|. Besides, an interpolating MLS is 
developed recently. For error analysis and applications to element-free Galerkin method 
see 


23, 


24|. 


The present work is based on the theory of Wendland and extends all the above results 
to a general case. All mathematical details are provided, special care is taken near the 
boundary, and lower bound for the minimum eigenvalue of the MLS local matrix is derived 
in general case, independent of the mesh-size. Besides, the analysis is presented for 
functions in fractional order Sobolev spaces. Finally an application to Galerkin methods 
for elliptic PDFs is investigated. 


2. MLS approximation 

Let O C for positive integer d, be a nonempty and bounded set. In the next 
section, more conditions on 17 will be imposed. Assume, 

X = {XI,X2, .. .,xn} C O, 

is a set containing N scattered points, called centers or data site. Distribution of points 
should be well enough to pave the way for analysis. 

Henceforth, we use for m € Nq = {n G Z, n > 0}, as the space of d-variable 
polynomials of degree at most m of dimension Q = A basis for this space is 

denoted by {pi,... ,pq} or As usual, B{x,r) stands for the ball of radius 

r centered at x. 

The MLS, as a meshless approximation method, provides an approximation Su,x of 
u in terms of values u{xj) at centers Xj by 

N 

u{x) fv Su,x{x) = ''^aj{x)u{xj), a; G 17, (2.1) 

where aj are MLS shape functions given by 

Q 

aj{x) = w(x,Xj) Xk{x)pk{xj), 

fe=i 
2 


( 2 . 2 ) 










where the influence of the centers is governed by weight function Wj(x) = w(x, Xj), which 
vanishes for arguments x, Xj € with ||x — Xj \\2 greater than a certain threshold, say S. 
Thus we can define 'Wj{x) = $((a: — Xj)/6) where $ : —>■ R is a nonnegative function 

with support in the unit ball B{0, 1). Coefficients Xk{x) are the unique solution of 
Q 

'^Xkix) ^ Wj{x)pkixj)pi{xj) = pe{x), 0 ^ < Q, (2.3) 

fc—1 j^J{x) 

where J{x) = {j : ||x — Xj \\2 ^ <5} is the family of indices of points in the support of w. 
In vector form 

a{x) = W{x)P^{PW{x)P^)~^p{x), 

where W{x) is the diagonal matrix carrying the weights Wj{x) on its diagonal, P is a 
Q X 4t^J{x) matrix of values pk{xj), j G J{x), 1 ^ k ^ Q, and p = {pi,... ,Pq)^ ■ In MLS 
one finds the best approximation to u at point x, out of with respect to a discrete 
H? norm induced by a moving inner product, where the corresponding weight function 
depends not only on points Xj but also on the evaluation point x in question. Note that 
A{x) = PW{x)P"^ is a symmetric positive definite matrix for all x G il. More details 
can be found in Chapter 4 of [l^. 

In what follows we will assume that <I> is nonnegative and continuous on and 
positive on the ball i?(0,1/2). In many application we can assume that 

4>(a:) = (/(||a;|| 2 ), a: G R"*, 

meaning that $ is a radial function. Here : [0, oo) —>■ R is positive on [0,1/2], supported 
in [0,1] and its even extension is nonnegative and continuous on R. 

If, further, (j) is sufficiently smooth, derivatives of u are usually approximated by 
derivatives of Su,x, 

N 

D°‘su,x{x) = ''^D°‘aj{x)u{xj), a; G H, (2.4) 

i=i 

and they are called standard derivatives. They are different from GMLS or diffuse deriva¬ 
tives |25| which are not the aim of this paper. 


3. Error estimation 

Since error estimates will be established using a variety of Sobolev spaces, we intro¬ 
duce them now. Let H C R'^ be a domain. For k G No, and p G [l,oo), we define the 
Sobolev space to consist of all u with distributional derivatives D°‘u G LPffX), 

jof] ^ k. The (semi-)norms associated with these spaces are defined as 

E > MHin) := f E 

\a\—k J \|ct|^fc 



\ i/p 
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The case p = oo is defined in the standard way 


sup ||T>“u||ioo(Q), ||u||vi/fc(o) :— sup ||T>“'u||ioo(n). 

|a|=fe |a|^A: 


For fractional order Sobolev spaces, we use the norms below. Let p G [l,oo), fc ^ 0, 
k G Z, and let 0 < s < 1. We define the fractional order Sobolev spaces to be 

the space of all u for which the norms below are finite. 


l^lwp*’+'’(0) ■“ 
ll'“llwp*^+'’(n) 


n Jn 


E 

Ja|=fe 

(ll“llwp‘’(n) 


\D°‘u{x) - D°‘u{y)\P 
|a; - y\^+P^ 

\ i/p 


i/p 


dxdy 


iw, 


p*^+'’(n)) 


The first step in deriving error estimates is to consider only local regions T) that are 
star-shaped with respect to a ball. A domain C is said to be star-shaped with 
respect to a ball B = B{y,p) = {x G : \\x — y\\ < p} if for every x G V, the closed 
convex hull of {x} U B is contained in 1). Let 


Pmax = sup{p : T) is star-shaped with respect to a ball of radius p}, 


then the chunkiness parameter of T) is defined by 7 = where d-p is the diameter of 
V. 

Approximating a function u G by averaged Taylor polynomials QmU G 


is discussed in [2^ Chapter 4]. The averaged Taylor polynomials are defined as follows. 
Let i? be a ball with respect to which V is star-shaped having radius p ^ ^Pmax- Then 

Qmu(x) := E A / - yT<^{y)dy, 

I ^ a! ds 

where p{y) > 0 is a C°° “bump” function supported in B satisfying both p{y)dy = 1 
and m ax^ ^ Cp~‘^. 

In [26| the Wp bounds on u — QmU are given for integer £ when u G and 

I ^ TO-1-1. A version of these results that applies when u belongs to 0 ^ s < 1, 

was proved in 271. An improvement of conditions (range of s) was discussed in [ 2 ^ by 
the same authors. 


Lemma 3.1. Let B be a ball in B such that B is star-shaped with respeet to B and 
such that its radius p ^ (l/2)pmax- Let QmU be the Taylor polynomial of order m of u 
averaged over B where u G Wf^~^^(B) for 0 ^ s < 1 and p G [l,oo). Let m > d/p for 
p > 1 and m ^ d for p = 1. Then there exists constant C = C{m, d,p, 7 ) such that 

||U - QmU||L~(D) < C (3.1) 


where dp is the diameter ofB. 
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We should note that the identity 


D°‘QmU = for all u € 


(3.2) 


which is found in 26|, Section 4], holds for |a| ^ m. Applying Lemma 13.11 on D'^u 
instead of u, using the identity (13.21) and the inequality |D“u|j^ic+ 3 -|o|p^ ^ \u\y^^k+s^^y 
we obtain 


Corollary 3.2. Letovs <1. For u G Wp+^{V), 

provided that m > \a\ F d/p for p > \ and to ^ |a| + d for p = \. 


(3.3) 


Analogous to the error bound (13.IL one can easily derive the bound for u — QmU 
when u G W™“'“®((D). We give the results in the following Lemma. 

Lemma 3.3. Let q G [1, oo], p G [1, oo) and a he a multi-index satisfying m > \a\F d/p 
for p > 1 and to ^ |a| + d for p = 1. With the notation and assumptions of Lemma \S. 11 
we have 


where C = C{m,d,p,q,a,^). 


Proof. Although the proof can be implicitly extracted from 27|, but we present it here 
for the reader’s conveniences. Let q G [l,oo). Using the definition of Sobolev norms, we 
have 


|| m - J^\D^iu-QmU)\'^dx 


l/3K|a| 


^ #{/3 G No : |/3| ^ |q;|} x vol{V)( max \\D^{u - Qmu)\\L^(^v)) 


^C{d,a)d^ max {u - Qmu)\\% 






< C(to, d,p, a, 7 ) rf^(-+-l“l+Ai/ 9 -i/p)) |„|9 


wF+‘’{v)- 


At the third line above, we use the facts that vo\{'D) ^ C^d^ and 

#{/3 £ : W\ < |o|) = i: (' +1:, ') = “) = OQaf). 

In the last line. Corollary 13.21 has been applied. Finally taking the g-th root of the both 
sides completes the proof with the new constant C = C{m, d,p, q, a, 7 ). The case q = 00 
can be proved by a similar argument (see also 1^, Proposition 11.29]). □ 
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Remark 3.4. In case s = 1, if we assume m + 1 > \a\ + d/p for p > 1 and m + 1 ^ |a| + d 
for p = 1 then estimations (EU and EH) are still valid, due to [^. The reader should 
be cautious that these error bounds can not be obtained by inserting s = 0 and replacing 
m by m + 1 in fractional cases, because the later produces Qm+iu. 


Up to this point, we reviewed some Sobolev error bounds for a function which is 
approximated by the averaged Taylor polynomial on a star-shaped domain. These bounds 
are usually used for analyzing the finite element method (FEM). Now we turn to the MLS, 
as a meshless approximation method, and employ the above bounds to analyze it. The 
final bound will be presented for functions in fractional Sobolev spaces. Although one can 
use the interpolation arguments (for example the “real” method based on K-functionals) 
to extend the integer order Sobolev spaces to fractional ones, here we follow the direct 
approach because all materials are provided via EH) and EH- 

First we introduce some other notations. For a set of points X = {xi,X 2 , ..., Xn} in 
a bounded domain C the fill distance is defined to be 

hx,Q. = sup min ||a; - Xj\\ 2 , 
xen KjAAf 

and the separation distance is defined by 

= \ Tam.\\xi-Xj\\2. 

A set X of data sites is said to be quasi-uniform with respect to a constant Cqu > 0 if 

qx < hx,n ^ Cquqx- (3.5) 

A set A = {ii,..., xx} C with N ^ Q is called P((j-unisolvent if the zero polynomial 
is the only polynomial from that vanishes on X. 

A set U C is said to satisfy an interior cone condition if there exist an angle 

0 G (0,7r/2) and a radius r > 0 such that for every x G fl a unit vector ^{x) exists such 
that the cone 


C'(x,^,6»,r) := {xty: y G W^, \\y \\2 = 1,2/'^C > cos6»,t G [0,r]} 


is contained in U. 

Assuming the compact set U satisfies an interior cone condition with radius r and 
ang le 9, and data site X C fl satisfies the quasi-uniform condition (13.511 . Wendland 
15l Chapter 4] proved that shape functions {aj(x)} from MLS approximation (12.11) 
provide a stable local polynomial reproduction of degree m on i.e. there exist constants 
ho. Cl, C 2 > 0 independent of X such that for every x G O 


1- = P{x)i for all p G 
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2- E^i l«i(2;)l ^ C'l 

3. aj(x) =0 if ||x - Xj\\2 > S = 2C2hx,a, 


for all X with hx,n ^ ^o- Constant C'l depends on the weight function </>, and constants 
C 2 and ho are 


C2 = 


16(1 + sin0)^TO^ 

3 sin^ 6 


ho = 


C2 


(3.6) 


He also proved that, if the weight function possesses k continuous derivatives then the 
approximant Su,x is also in C^. Using the above properties, he proved the error bound 



where |u|c'"+i(f2*) '= in which Vt* = \J^(iQ,B{x,C 2 ho) can be 

obviously larger than the exact domain U. Here is the classical derivative operator 
on space The results of the present paper (in a special case) extend this bound 

for functions u in fractional order Sobolev space W^“''®(H). 

Taking derivatives of order a, |a| ^ m, under some mild conditions, we can show that 
functions {Il“aj(a;)} in approximation (12.41) form another local polynomial reproduction 
in the following sense: there exist constants hg, C '2 > 0 independent of X such that 
for every x G U 

1- D'"o,j{x)p{xj) = D°‘p{x), for all p £ 

2. l^“«i(^)l ^ Ci,„h-'“', 

3. D°‘aj{x) =0 if ||x - Xj||2 > ^ = 2C2hx,n, 

for all X with hx.n ^ hg. The first and the last items are immediately followed from 
the previous local polynomial reproduction system. But proving item 2 invites more 
challenges. First we prove the following straightforward result. 

Lemma 3.5. Let X = {xi,... ,xn} C U has fill distance hx,n- Suppose that function 
4 > : [0,oo) —>■ R, is supported in [0,1] and its even extension belongs to C’"(R) form £ Ng. 
Then for Wj{x) = <p{\\x — Xj\\2/5), x £ and |a| ^ m we have 

|L»“u;,(x)KC„/i-|“', Vx£H, j = l,...,N, (3.7) 


provided that 6 = 2 C 2 hx,n- 


Proof. Since ()) is a compactly supported and C"" function, derivatives of (j) up to order 
m are continuous and bounded. The absolute value of D°‘wj, has a bound with a factor 
times derivatives of 4>. This immediately gives the desired bound for sufficiently 
small hx,n- □ 


7 




The MLS approximation can be implemented in a more stable fashion, if a shifted 
and scaled polynomial basis function is used as a basis for In this case, we use the 
basis 

f (x- z)“ 1 


- I“l 


(3.8) 




where z is fixed and depends on the evaluation point to be considered. If x is the 
evaluation point, the best result will be obtained if we finally set z = x. In fact, MLS 
uses different bases for each evaluation point. We can do this, because the formulation of 
MLS approximation and equations (12.21) and (12.31) are independent of the choice of basis 
functions. Thus the MLS shape functions can be written as 


aj{x)=Wj{x) ^ K{x) 


{Xj - zf 
Ja| 


j = l,2, 


,iV, 


(3.9) 


where \a (x) is obtained by solving the positive definite system 


A(x)A(x) = p(x), Pa{x) = 


(x - z)^ 


^l«l ( 

^5^ ^ 0^|a|^m 


where 


N 


Ai,kix) = y^Wj{x)pi 


J=1 




Pk 




i,k = 1,... ,Q. 


(3.10) 


(3.11) 


N 


Since ij) is supported in the unit ball, we used the summation index instead of 

j=i jeJix) 

in the above formulation. Since the set point X satisfies the quasi uniform condition 
(1531) . the number ^J(x) of points in J(x) can be bounded independent of hx,(i |15| . In 
fact, for Xi,Xfc € B(x,S) and Xi ^ Xk the balls B(xi,gx) and B{xk,qx) are disjoint. All 
of these balls with Xj € J{x) are contained in the ball B{x, qx + <5). It is clear that 


vol (u B{xj,qx)'j i^vo\{B{x,qx + S)), 


jeJix) 

which simply gives 4fJ{x)qx ^ (<5 + 9x)^- Using the quasi-uniform condition and (5 = 
2 C 2 hx,n, we have 

#J(x) ^ (l + 2 C 2 Cqu)''=:U#. 

Lemma 3.6. If the weight function (f satisfies the assumptions of Lemma \3.5[ then for 
a fixed but arbitrary evaluation point x G H we have 

\D°‘A{x)\ ^ Va with |a| ^ m, 

where Ca is a constant matrix independent of hx,n- 
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(3.12) 








Proof. Equation (jS.lip gives 


N 


D'^Ai^kix) = '^D°‘wj{x)p^ 
i=i 




Xn — X \ I Xi — X 

Pk 




Evaluating at x, taking absolute value from both sides and using ||a:j —iE ||2 ^ <5 = 2C2hx,n 
we obtain 


N 


\D^AA^)\ < \D^Wj{x)\ ^ 

i=i 


This completes the proof. 


□ 


Since A{x) is positive definite for all x G il, all eigenvalues are real and positive. If 
the basis (Ira is employed, we can prove that the smallest eigenvalue of A{x) has a lower 
bound away from zero and independent of hx,Q.- Proving this assertion helps us to find 
a bound for \D°'A~^{x)\. First, recall 

Ainin(^(a;)) = min , (3.13) 

i;eK«\{0} V 


for symmetric matrix A{x). Since A(x) is also positive definite, we necessarily have 
Amin(^(a;)) > 0. To bound Amin we follow some parts of Melenk’s argument [21| and the 
concept of norming sets presented in [Appendix .![ 


Lemma 3.7. Suppose that the bounded set il C satisfies an interior cone condition 
with radius r and angle 6 G (0,7r/2). Let X = {cci,..., xjv} C 17 has fill distance 
h = hx.Q, and satisfies h ^ r/C 2 ='■ hg where C 2 is defined in (13.611 . Suppose that 
d = 2 ( 72/1 is the size of supports of the weight functions, and for a fixed but arbitrary 
X in LI the set {xj G X : j G J(a?)} is Wf^-unisolvent. Then there exists constant C\ 
independent of h such that 

Amin(^(a;)) > (7 a > 0, 

provided that the shifted scaled basis functions dSH) are employed. 


Proof. Let v* G v* ^ 0, be a vector at which the minimum in (13.1311 occurs for 
X = X. Define 

1 ( 3:1 

Now using f|3.11p we simply have 


N 

V*^A{x)v* = ^Wj{x)TT^{Xj). 
1=1 
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Since {xj : j G J{x)} is PJ^-unisolvent, the functionals Z = {Sx^ '■ j G J{x)} form a 
norming set for i.e. there exists an injective mapping T : P^ —>• T(P5^) C 
where T{p) = {p{xj ))(See [Appendix . l] ) . Set Z allows us to equip PJ^ with an 
equivalent norm via the operator T. We define the norm on by 

WpWIw ■= 

j€J(x) 

and the norm on PJ^ by the infinity norm. Using the properties of norming sets and 
setting p = TT, we have 

^ ||J1—1|| ll’^lloO,B($‘,(5) ■ 

Obviously, the set dB(x, h) = {x : ||a; — ai|| = h} is a subset of B{x, S). Definition of tt in 
(13.1411 ensures that the values of tt on dB{x, h) are independent of h, because h will be 
canceled from the numerators and the denominators. In fact ||7r||oo,aB(s,;i) is bounded 
from below by a positive factor times ||v*||i := X]|a|^m have 

l|7r||oo.B($.<5) ^ |K||oo.aB(£,/i) ^ C' 7 r||u*||i > 0^11^^*112, 


where O^ is the mentioned factor which is independent of h. The last inequality follows 
from the standard relations between one and two norms in R*^. 

It remains to bound ||T“^||. By assumptions, ft satisfies a cone condition with angle 
9 and radius r, and h < r/C 2 - The later gives 6/2 ^ r. Of course the cone condition 
will be obeyed if we use any radius less than r. Thus for every x G D, there exists a cone 
C{x) = C{x,^,9,S/2) C V, (1 B{x,S/2), and using Lemma p\ppendix .2| there exists a 
closed ball 

B = Bix,pS)cCix), 

We are going to prove 


lblloo,B(5r.5) < C'llbll^ B ^ C'b(a;fc)|, Xfc G 5 c B{x,6/2), 


(3.15) 


for all p G P^. Using Lemma Appendix .4[ the first inequality satisfies with Ci = (|)™- 
To prove the second, let Ibiloo b ~ \p{xm)\, for xm G B. Since ball B itself satisfies 
an interior cone condition with radius p5 and angle tt/S, Theorem [Appendix . 3| can be 
applied provided that h ^ ^ . One can easily check that this condition is always 


satisfied because 5 = 2 C 2 h = ^ 

^ 3sin^ 0 


4(l + \/3/2)m2 ■ 

-l^2 2 

. The proof of the mentioned Theorem shows 
that there exists a point Xk & X (1 B such that \p{xk)\ > \\p{xm)\ = ^Ibiloo S' Thus in 
(13.1511 the constant C can be chosen as Ce,m '■= )™. Letting 


uimin := min ^(s), 

sG [0.1/2] 
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we can write 


N 


W^min|b||L,B(£,5) < C!lmWk\p{xk)\^ ^ ^ ) P = C'^.mlblll.^ 

i=i 


which immediately gives 

Summarizing all, we have 






X V*A{x)v* ||7r||i,^ WminC^ ^ 

Ai„in(^(a:)) = - TtV-;— = H ^1,9 ^ n^.nn ^ ^yj- =: ^A, 


V*T'v* 


IMIi ^ \\T-^"^ 


ri2 

^e,m 


which completes the proof. 


□ 


Remark 3.8. The unisolvency condition in Lemma 13.71 is a mild condition, because in 
most cases if satisfies a cone condition then fl fl B{x,6), x G il, also satisfies another 
cone condition and for sufficiently small hxnB,nnB Theorem Appendix TS ensures the 
unisolvency. 

Remark 3.9. The role of “shifted” and “scaled” basis functions (13.81) is crucial to bound 
Amin away from zero and independent of hx,n- Otherwise, experiments show that Amin 
tends to zero when hx,n 0. See section 6 of for numerical results. 

Lemma 3.10. With the notation and assumptions of Lemmas and \3 . 51 we have 

\D°‘A~^{x)\ ^ Cah~'^°‘\ Va with |Qf| ^ m, (3.16) 


where Ca is a constant matrix independent of h = hx,n- 

Proof. Taking the derivative of both sides of the well known relation A~^{x)A{x) = / 
and evaluating at x, we obtain 

D^A-\x) = -A-\x) {^\d^-^A{x)D^A-\x). (3.17) 

Pita 

Induction on |/?| can be used to prove the desired result. The first step in induction is 
considered by choosing (3 = ej, the unit vector with 1 in j-th place. We simply have 
D^^A~^ = — A~^A)A~^. Since \ej\ = 1, equation (13.121) yields 


where Ci is a constant matrix. From matrix computations, there exists a constant matrix 
C G R'SxQ |A“^| ^ C'||A“^|l 2 holds. Since A is symmetric positive definite, we 

have 




Ch-^ C 

AST 
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which yields the starting point for induction. Now we suppose that \D^A ^ Ch 1^1 
holds for all /3 ^ a and /3 ^ a. Employing (I3.12|) . equation (13.171) gives 


\D^A-^\ ^ |A-i| ^ ^ 

/3^a 

/3/a 

which completes the proof. □ 


Theorem 3.11. The shape functions aj, j = from the MLS approximation 

possess the following stability condition for |a| ^ m, 

N 

VxGfi, 

i=i 

where Ci^a is independent of data site X, provided that all assumptions of Lemmas \3. 7| 
and \3.5\ are satisfied. 


Proof. Let h = hx.Q. First, from (13.1011 we have A(a;) = A ^{x)b{x) and thus for a fixed 
but arbitrary x G LI 

D°‘X{x) = ^ ('^'\d°'~'^A~^(x)D'^p{x), Ma with |a| < to. 

Since z = x, obviously all entries of vector D^p{x) are zero except the ry-entry which is 
i.e. D^p{x) = where e,, is a unit vector with 1 in 77 -th place. Now 

using (13.161) we can write for a constant matrix Ca,r) 



where the vector Ca is a bound for 
sides of equation dsai one obtains 


Ca^nS-ri- Now, taking the derivatives of both 


iA“a, 



|/3|^m 



Evaluating at x, applying the bounds of D'^\p{x) and D°‘ ‘^Wjfx) and using the fact 
that \xj — x\^ < hl^l, we finally have 


7]^a I |/3|^m 


rj^a 

if 


□ 


which completes the proof. 
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Theorem 13.111 establishes the second property of the local polynomial reproduction 
system {D°^aj}. This will help us to estimate the error function in MLS approximation. 
First we note that a region with a Lipschitz boundary automatically satisfies an interior 
cone condition. More details can be found in 


29|. 


Theorem 3.12. Suppose that 17 C is a bounded set with a Lipschitz boundary. Let 
m be a positive integer, 0 ^ s < 1, p € [l,oo), q € [l,oo] and let a be a multi-index 
satisfying m > la] + d/p for p > 1 and m ^ |q;| + d for p = 1. If u € (LI), 

there exist constants C > 0 and ho > 0 such that for all X = {xi,... ,a:Ar} C 17 with 
hx,n ^ min{/io,l} which are quasi-uniform with the same Cqu in \3.5\) . the estimate 


u — Su,X 






ws 


Tf2)> 


(3.18) 


holds. Here (a;)+ = maxjx, 0} and Su,x is the MLS approximation of u on data site X 
in which the corresponding weight function satisfies the assumptions of Lemma \3.5[ and 
the shifted scaled basis polynomials (13.81) are employed. 


Proof. Since 17 is bounded and has a Lipschitz boundary, we can use the continuous 
extension operator 


En : LF™+*(17) ^ 1 < p < oo, 

to extend any u G LF™+®(17) to a function v := Equ G with u|n = u. Since 

the extension is continuous, we have 


||u||^m+3(jjd) ^ C||u||^y™+*(Q). (3.19) 

The case s = 0 was constructed by Stein and works also for p = oo. DeVore and 
Sharpley have proved this extension for the fractional order spaces. 

First we prove (13.1811 for q G [1, oo). The case q = oo will be discussed later. Let the 
Lipschitz domain 17 satisfies a cone condition with angle 9 and radius r. Assuming ho = 
r/C-i in (13.61) . we first bound the error over subdomains Bk = B{xk, d) (~|17, k = 1,..., iV, 
for 5 = 2(72 hx.n where hx,n ^ min{/io, 1 }. At the end, we will extend the error bound 
over entire 17. Let Vk = B{xk,‘25), k = 1,...,N. Clearly, Vk ^ 17 in general. But Vk 
is star-shaped with respect to a ball B C T>k with chunkiness parameter 7 = 2 . Now 
let p = QmV G be the Taylor polynomial of degree m of u on Vk averaged over B. 
The reader should care about the letter p, which has been employed for both polynomial 
and Sobolev notations. Using the properties of the stable local polynomial reproduction 
{ttj}, we can write for x G Bk 

N 

u{x) - Su,x{x) = u{x) -p{x) +'^aj{x){p{xj) - u{xj)), 

t=i 
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and in w}^'^ norm, 


N 


||M-s«.x||^i<.|(g^) < \\u - + \\J2aj{-){p{xj) - u{xj)) 


i=i 




(3.20) 


Using the facts that Xj € T>k and z;|q = u, the second norm on the right-hand side can 
be bounded as below 


N 


N 


i=i 




l/3|^|a| i=l 


Q 

L^{Bk) 


N 


|/3|^|a| "" j=l 


l/3K|a| 


^ C'c4,Ua'^lk-p|li 


where we use 

N 


Y1 

i=i 


N 


L«>(X)fc)’ 


< / f W |D^ai(a;)| |p(a;,) — u(a;,)|') dx 

i’(B.) ^ II ^ ^ ly 

Af 9 /■ 

< lk-p|II=o(i,,,)(max^ |D'^aj(a;)|) dx, 


■i=i 


which together with Jg^ da; = vol(i3fc) ^ cdg^, gives the inequality in the second line. 
The inequality in the third line follows from Theorem 13.111 The last estimate satisfies 
because hx.n ^ 1 and |a| ^ |/3|. Thus from (13.201) we can write 


\W- Su,x\\w!,‘^'(b^) ^ lk-)PllH^i“i(B,) +'^«dg///i^|“'||z;-p||icx=(i,,). 

To bound the both terms on the right-hand side of inequality above, first by (EH) we 
have 

Ik-P||L~(I5A,) < 

Then, since Bk CVk, (EH) leads to 

\\^-P\\wl-\B,)^\\^-P\\wi-'iV,) 

,mS 


^ ,m+s-\a\+d(l/q-l/p)t 

^ O C Ti/m+s/ 


m^^+^CDk)- 

If we assemble everything up to this point and use the facts that dg^ ^ 25 and d-p^. = 4(5 
we get 

Ik - llvvi<»i(B,) ^ (3.21) 
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Now we should extend this bound over entire Since 6 = 2C2hx,Q. and C 2 ^ 1/2, 

for every x £ there is a center Xj G B{x,S) fl 11. This clearly shows H C 

=: fl*. First, since Vk C fl* we have 


N N 

X! - - J 

k—1 k—1 \^\—m ^ 




EE 


N 

E 

k—1 \^\—m 


T>k 


f \Df^u{x) - D^u{y)\P 

r \Dl^u{x) - Dl^u{y)\P 
/n- k - y\<^+P^ 


dxdy 


dxdy 


N 


= E / 


\k^l 


f \D^u{x) — D^u{y)\P 
/o. \x - y|^+P« 


dxdy, 


where xv denotes the characteristic function of the set V. Note that n{x) := J2k=i XVk {^) 
is the number of subdomains Dk containing x. This function can be bounded by a con¬ 
stant because X is a quasi-uniform set. In fact n{x) is the number of points Xk located 
in the ball B{x, 26). Since this ball is contained in a cube of side-length ASjVd, we can 
write 


n{x) ^ 


/ 4(5 


\Vdq:^ 


€ 


8C'2hx,nCqu 

Vdhx.Q, , 


^Cq\lC'2 


Thus we have 


N 

El 

fc=i 


'w; 


“(•DO 




^E 


= C\-'P 


\P\=m 

v\ 


r r \DPu{x) - DPu{y)\P 

/o- in* 


dxdy 


'w^+‘{n*y 

Now applying p.2ip and the above bound we can write 

/ N \ 1/9 


\fc=i 




N 


1/9 


Ei^^i 




Wp”*+"(I5fc) 
/ N 


1 /p 


€ 


(7/j™+®“l“l+‘^(i/9“i/p)_/y(i/9-i/p)+ j 1^1 


wy 


Kk=l 


“(DO 


^ C/ i^x.Q n TP o 


x,a 


m; 


“(f2*) 




^ r'1„| 

^ '^x.n \^\w, 


lVp'"+Tf2*) 
+“(R'i) 


w; 
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The bound on the third line above follows from standard inequalities relating p and q 
norms on finite dimensional spaces where (a;)+ = max{a;, 0}. In the fourth line, to bound 
N by the fill distance, let da be the diameter of ft. Since ft is bounded, there exists a 
cube of side length da/Vd that contains 17. Thus 


^ydqx 


^ f c^udn 


^ i/dh 


x,n' 


x,n- 


In the fifth line, we have used the identity d(l/g—1/p) —d(l/g—l/p)_|_ = —d(l/p—l/q')^. 
Finally, we invoke the norm equivalence property (13.191) to get the final bound 

The case q = oo can be proved in a similar way, because dSH) can be used for g = oo to 
bound the first term in (13.201) . and the second term can be simply bounded by 


N 


|^aj(-)(p(xj)-u(xj)) 

t=i 

The reader can continue the proof to get 




wls'iBk) 


II ^ im+s—\ot\ — d/p 

u — ^ ^ ^X.Q 


Iw”+Tf2)- 


□ 


Remark 3.13. According to Remark l3.41 one can easily proceed with the proof of Theorem 
13.121 (by doing some modifications) to get the estimation 


«(3.22) 

provided that m + 1 > |a| + d/p for p > 1 and m + 1 ^ |a| + d for p = 1. 




4. Application to Galerkin method for PDEs 

As an application, we consider the second order elliptic partial differential equation 

“ ^ ^ +c(x)M(a:) = f{x), x€ft, (4.1) 

Ki,(x)--— (x)ni(x) + b(x)u(x) = g(x), x € dft, (4.2) 

axi 

i,j=i ■> 

where 17 is a bounded domain with Lipschitz boundary 517, and Kij,c € L°°{ft), / G 
L^(17), aij,b G L°°{dft), g G L^(dft) and n is the unit normal vector to the boundary 
517. Matrix K(x) = («)y(x)) is assumed to be uniformly elliptic in 17, i.e. there exists a 
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constant 7 such that for all x € and all a G we have a^K{x)a ^ TllalU- Moreover, 
we assume c ^ 0 and 6 ^ 0 , and at least one of them is uniformly bounded away from 
zero on a subset of nonzero measure on Q or dfl, respectively. Under these assumptions 
the weak form of equation (14.11) together with boundary condition (14.21) is a{u, v) = £{v) 
where a{u, v) : x Wl{^) —>■ R is a coercive and continuous bilinear form defined 

by 


a(M, v) 



du dv 

—"5 - 

oxj axi 


dfl ' 


buv dr, 


Ida 


and £ : W 2 (fl) —>■ R is a continuous linear functional defined by 


£{v) = j fvdQ+ j gvdr. 

Jn Jan 

Using the Lax-Milgram theory, the corresponding variational problem 


find u G Wl{£X) such that a{u,v) = £(v), for all v G Wl(Vt), (4.3) 


admits a unique solution u and the solution is continuously depended on data £. This 
problem has been analyzed in using radial basis functions interpolation. 

To find the numerical solution we use the same Galerkin method as in the classical 
finite element method. The approximation solution is sought in a subspace generated by 
MLS shape functions. We define for quasi-uniform set X = {xi ,... ,xn} C U 

Vn =span{ai,a 2 ,...,aAr} 
as a subspace of Wl{^) and solve the discretized problem 

find un G Vn such that a{uN,v) = £{v), for all v G Vn- (4.4) 


Of course this step concerns the computation of domain and boundary integrals, which is 
the most difhcult stage of the procedure. But we assume that all integrals are computed 
accurately and seek a bound for the error ||m — UAr||wi(n) for fbe function u G 
where m > 1 -I- d/2 and 0 ^ s < 1. Our analysis allows to consider functions that are 
less smooth than the functions in First, recalling the Cea’s Lemma we have 

llw-Mivllwi(n) ^ C* mf ||u-^^||wi(n): 
where C is a generic constant. Since Su,x G Vn, we obtain 

\\u — UN\\w^{n) ^ C” ||m — Su,x\\w}{n), 
which leads to the following corollary. 
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Corollary 4.1. Let Lt C be a hounded domain with Lipschitz boundary, and m be 
an integer satisfying m > 1 + d/2 and let s G [0,1). Suppose that u G is 

the solution to the variational problem (USD and un G Vn is the solution of discretized 
problem (HU) where Vn is constructed by the guasi-uniform set X = {xi,... ,Xn} C LI, 
the weight function (j) satisfying assumptions of Lemma \S.5[ and the basis functions (IH.8I1 . 
Then there exist constants C and ho such that for all set X with hx,n ^ min{/io, 1} the 
estimation 

||u — UAr||^v^i(f2) ^Ch^'^ 111114 ^ 2 ”“+® (n) 

holds. 

Finally according to (I3.22|) and discussions before Corollary 14.11 if integer m satisfies 
m > d/2 and u G then the error bound 

||u— UAr||vy2i(Q) < 

holds. 

The orders are the same as those for classical finite elements. In both cases we can 
use the technique of Nitsche to estimate the error in L^-norm. 

5. Numerical examples 

Since there are extensive numerical examples in literature, here we will restrict our¬ 
selves to a couple of examples, in which we will concentrate on the predicted orders of 
the errors in (13.1811 and (13.2211 . 

We consider the following example 

u{x) = ||a;|| 2 , X G Tl C 

where A is a real parameter and is a bounded region around the origin. It is well known 
that 

u e Wp (O) A > T — d/p. 

We let p = 2 and If = [—0.5,0.5]^ C and we assign two values 1.5 and 3 to A. 
According to the theory, in the first case we set m = 2 and examine (j3.18|l . and in the 
second case we set m = 3 and examine (13.2211 . In both cases a regular mesh distribution 
with the fill distance h is used as a set of centers, the compactly supported C'^ Wendland’s 
function (j){r) = (1 — r)[j_(35r^ -|- I8r -|- 3) is employed as a weight function, and 5 = 2mh 
is used as a support-size. Results are presented in Tables [T] and [5] for q = 2, 00 , and 
different order derivatives a. The L 2 -errors are computed using a (200 x 200)-point 
Gauss-Legendre quadrature, and Loo-errors are computed on a very fine regular mesh of 
size hs = 0.005. 

As we can see, the experimental results confirm the theoretical bounds. 
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Table 1: 

Orders for A = 

1.5 and m = 2 



L 2 


^00 


h 

a = (0, 0) 

a = (1,0) 

a = (0,0) 

a = (1,0) 

0.1 

0.05 

2.56 

1.51 

1.50 

0.50 

0.025 

2.52 

1.50 

1.50 

0.50 

0.0125 

2.56 

1.49 

1.50 

0.62 

Theory 

2.5 

1.5 

1.5 

0.5 


Table 2: Orders for A = 3 and m = 3 



L 2 



LqO 



h 

a = (0, 0) 

a = (1,0) 

« = (2, 0) 

a = (0, 0) 

a = (1,0) 

a = (2,0) 

0.1 

0.05 

3.75 

3.08 

2.06 

3.00 

2.03 

1.00 

0.025 

3.86 

3.01 

1.99 

3.00 

2.00 

1.00 

0.0125 

3.88 

3.00 

1.93 

3.00 

2.06 

1.00 

Theory 

4 

3 

2 

3 

2 

1 


6. Appendix 

We restate the following definition, lemmas and theorem from Chapter 3 of the book 

0 . 

Definition Appendix .1. Let F be a finite dimensional vector space with norm || • ||y 
and let Z ^ V* (the dual space of V) be a finite set consisting of N functionals. We 
will say that Z is a norming set for V if the mapping T : V ^ T{V) C defined by 
T{v) = (z(v))zez is injective. T is called the sampling operator. 

If Z is a norming set for V, then T~^ : T(V) —>■ V exists and we can simply show 
that 

llTulk. ^lirillHlv, |H|v^||r-i||||rn||«., 
which means that || • ||v and ||T(-)||RAr are equivalent norms on V. 

Lemma Appendix .2. Suppose that C = C{x,^,0,r) is a cone. Then for every 
0 < h ^ r/{l + sin0) the closed ball B = B(jj, hsinO) with center y = x + hf and radius 
hsinO is contained in C{x,£^,9,r). 

Theorem Appendix .3. Suppose that 12 C is compact and satisfies an interior 
cone condition with radius r > 0 and angle 6 S (0,7r/2). Let m G N be fixed. Suppose 
h > 0 and the set X = {xi,X 2 ,... ,xn} C 12 satisfy 
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( 1 ) 


r sin 9 

4 (l+sin 9 )'^m^ ’ 


(2) for every B{x, h) C ft there is a center xj G X H B{x, h); 


then Z = {i5xi,..., is a normig set for ^^6 inverse of associated sampling 

operator is bounded by 2. In fact for every p G there exists Xk G H D X such that 
\p{xk)\ ^ ^|b||oo,n- If h = hx,n, the second item is automatically satisfied. 

Note that, the functionals Z = {5xi, ■ ■ ■ form a normig set for if and only 

if X is ¥l^-unisolvent. 


The following Bernstein inequality can be easily proved by using the one dimensional 
Bernstein inequality 


||p||oO,( —p,p) ||p||oO,(—1,1) ! y P G 


ll 

m ■ 


Details of the proof can be found in 2l|, Lemma B.4]. 

Lemma Appendix .4. Assume that Bi and B 2 are two balls of radius pi and p 2 , 
respectively, and Bi C B 2 c¥‘^. Then 


||p||oo,B2 ^ 


2p2 

Pi 


Iblloo.Bi, Vpe 


■jid 

■ m ■ 
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